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Abstract 



We discuss how the proton structure function F2(x,Q 2 ) is described in the HERA 
kinematic range by double asymptotic expressions for low x and large Q 2 . From a NLO 
double asymptotic approximation of recent data from the HI experiment at HERA we 
extract the strong coupling constant a s (M|) = 0.113 ± 0.002(stat) ± 0.007(syst). The 
additional theoretical error can be as large as 0.007. 



Previous fixed target lepton-nucleon scattering experiments [jl], |2| measured the proton structure 
function F2(x, Q 2 ) at Bjorken-x > 0.01 and were sensitive to the valence quark content of the nucleon. 
HERA experiments collect data at values of x low enough to neglect valence quarks and at values of 
Q 2 large enough to apply perturbation theory. Here Q 2 is the four-momentum transfer squared from 
the incoming lepton to the proton. The proton structure function Fi strongly rises with decreasing x 
over the whole kinematic range accessible to the HERA experiments HI and ZEUS ||. Such a rise 
was predicted more than twenty years ago || from the leading order renormalization group equations 
of perturbative QCD. Ball and Forte pursued these ideas || and proposed a way to demonstrate that 
the low x data at HERA exhibit scaling properties dominantly generated by QCD radiation. They 
gave an expression for F2(x,Q 2 ) in the double asymptotic limit of low x and large Q 2 . The recent 
F2(x, Q 2 ) measurement of HI is in accord with such scaling behaviour. Hence these data are expected 
to be sensitive to the fundamental QCD evolution dynamics, and not to depend on unknown starting 
distributions. The formalism of the double asymptotic analysis was extended to two loops in 

We define the logarithmic QCD evolution variables 

t = \n(Q 2 /K 2 ) and £ = ln(l/z) 
where A is the QCD mass scale, and the evolution length T of as(Q 2 ) from a starting point Qq to Q 2 

T = In (a s (Q 2 )/as(Q 2 )). 

To leading order T is simply T = 4ir/bo J* t * dt a>s(i) = ln(i/io)- The T dependence of F2 allows to 
determine as- Here bo = 11 — 2nj/3 is the leading coefficient of the renormalization group equations, 
with rif the number of flavours. 

For large t and £ and an Fi{x, Q 2 ) which at Qq is not too singular in x the NLO double asymptotic 
expression for F2(x,Q 2 ) is Q 

F2 ~ Np (1-fNLo) e X p(2 7v /^-<5T+ilnT-^m£). (1) 

with 

7 2 = 12/6o 5= (ll + 2n / /27)/6 - 

The normalization coefficient is Np = 5nj/324. For nj = 4 this gives 7 = 6/5 , 5 = 4/3 and 

Np = 0.038. The NLO correction term Jnlo is 

Inlo = ^^K«s(Qo) " «s(Q 2 )) - 13«s(Q 2 )] 

with e = (206n//27 + 66i/&o)/&o an d b\ = 102 — 38nj/3. The leading order formula is recovered 
setting /nlo = ^ an d b\ = 0. Defining the leading exponent as A£ = I^^JTT and the subleading term 
as a = —5T + l/41nT — 3/41n£ we rewrite F2 as 

F 2 ~ x~ x e a . (2) 

The leading term in the double asymptotic formula for F2 corresponds to the double leading log 
approximation DLL Q of the DGLAP equations Q. It generates the growth of the structure function 
with falling x proportional to x~ x . The subleading term falls with x but slower than the leading term 
grows. 

Expression (1) for the structure function is a prediction of QCD in the limit of small x and large 
Q 2 . It does not depend on the shape of the parton distribution functions at Qq if these are sufficiently 
soft. The starting point Qq of the QCD evolution is a free parameter. Conventional QCD analyses of 



structure functions [|10| , |l lL 12, |13[ fit input parton distributions at some Qq and evolve them to larger 



Q assuming a singular behaviour F2 ~ x of the input distributions. On the contrary, the approach 



of GRV [14] starts with non-singular, valence- like input distributions at some low Qq and generates 
the steep rise of the F 2 from QCD dynamics. 



Independently of a set of input parton distributions we now determine A and as fitting expression 
(1) for F 2 (x,Q 2 ) to the latest measurement of the proton structure function by the HI experiment 
JH at HERA. When passing the heavy flavour thresholds we adjust A according to [15]. In order to 



reach the asymptotic region in T and £ we demand Q 2 > Q 2 m i n = 5 GeV 2 > Q\ and x < 0.1 which 
corresponds to £T > 2 approximately. The fit using the total errors gives the central value A = 248 
MeV with a x 2 /ndf = 122/145 and Qq = 1.12 GeV 2 . The A symbol is taken as an abbreviation for 
A^L.[[lll 16|. The result is shown as a full line in Fig. 1 and describes the data for Q 2 > 5 GeV 2 
reasonably well. A fit of the LO expression of equation (1) to the F 2 data is found to drop too fast 
with x. The x 2 /ndf increases by a factor of two. 

The statistical error is 5A stat = 23 MeV from a fit with the statistical errors only. Using the full 
systematic error matrix of HI we then vary each of the measured quantities entering the F 2 analysis. 
From the variations of the A determination with the varied F 2 we get the full systematic error 8A sys t = 
83 MeV . In order to estimate the systematic error coming from the kinematic cuts we also varied 
Qmin f rom 5 to 60 GeV 2 and asked for x < 0.05 or x > 0.001. We also replaced the upper x limit by 
demanding that £T > 2. The resulting A variation is 46 MeV. The A determination translates into 
a 5 (M|) = 0.113 ± 0.002(stat) ± 0.007(syst). 

The heavy flavour treatment causes a theoretical uncertainty. Since we work at Q 2 above the charm 
threshold only the bottom quark mass m& enters the uncertainty. We chose an effective = 4.74 



GeV according to [17|. Taking mj, = 4.50 GeV instead changes A by 13 MeV. Using a different 
prescription 18] for the transition between AW' and AW -1 ) changes A by 7 MeV. Changing as 
an extreme assumption the mass scale from Q 2 to W 2 changes A by 15 MeV. The heavy flavour 
treatment thus yields a systematic error on as of 0.002. In order to check the influence of the 
asymptotic expression (1) on the as determination we also fitted it to the HI QCD fit || and to 
a set of GRV parametrizations with varying A [19| for x > 0.001. For the Qq obtained above the 



A values from these fits follow the variations of the values used in the QCD evolutions of these parton 
densities and differ by at most 34 MeV from the input values. This corresponds to an as error of 
0.003. The largest theoretical error, however, is anticipated to come from higher order perturbative 
corrections and subleading terms. This error is quoted to be in the range of 0.004 |2^] to 0.007 | [21f , 
We conservatively quote our result with a theoretical error of 0.007. 

We finally determine a QCD inspired parametrization of F 2 which is valid down to the lowest 
Q 2 values. Since the NLO fit shows that the structure function has not evolved steep enough below 
Q 2 = 5 GeV 2 we attempted to fit also the constant of the leading exponent in equation (1) fixing 
F 2 = ax~ with A = c-y/T/£. Including now the data at 1.5 < Q 2 < 5 GeV 2 we get values for a and 

c very close to Nf and 7 for four flavours. Hence we fix a = Nf and c = 7 and fit F 2 = NfX~^^/ T ^ 
with rif = 4 for x < 0.1. The result is shown in Fig. 1 as a dashed line. The values of the only two 
fit parameters are Ql = 0.365 ± 0.026(stat) ± 0.048(syst) and A = (243 ± 13 ± 23) MeV. The two 
parameters resemble the GRV Qq scale and the QCD parameter A . The x 2 / n df of the fit using full 
errors is 109/167. Leaving the flavour number rif free gives very similar Qq and A values, yielding 
rif = 3.83 ± 0.18 ± 0.34 with almost the same x 2 / n df = 108/166. For x < 0.1 this fit describes the 
F 2 data surprisingly well with the two parameters Qq and A only. 



Summary 

For Q 2 > 5 GeV 2 a NLO double asymptotic expression for F 2 (x,Q 2 ) describes the HERA data well. 
We use its Q 2 dependence to extract as{M\) = 0.113 ± 0.002(stat) ± 0.007(syst). We quote a theory 
error of 0.007 which is subject to uncertainties. We found a simple QCD inspired parametrization 



F 2 = Vje _7 v / ^V€ for the proton structure function which is valid over more than three orders of 
magnitude in x and Q 2 including the region Q 2 < 5 GeV 2 and depends on two fitted scales Qq and A 
only. 
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Figure 1: The proton structure function F2(x,Q 2 ) as measured by the HI experiment at HERA 
together with a fit to the NLO double asymptotic expression (1) (full line) for Q 2 > 5 GeV 2 and with 
a fit to the modihed DLL expression F 2 = Nfe^Vr/l (dashed line) in the full Q 2 range. 



